Short-range defects contribution to the monolayer graphene resistivity 
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The derived by us earlier electron scattering matrix for the short-range defects in monolayer 
graphene is applyed to description of the resistivity electron density dependence. It is argued that 
large charged defect density is unlikely in the suspended graphene and graphene lying above the 
trench. 

I. INTRODUCTION 

Electronic states in monolayer graphene with point defects were considered in 0], Q and in our works Ql, 
0. The main novel element in our works is the band asymmetry of the defect potential in the Dirac equation. 
This asymmetry appears naturally if the defect violates the symmetry between the sublattices. This asymmetry 
was described in terms of equivalent superposition the band-symmetric potential and the mass (gap) perturbations. 
Characteristic equations for the energy eigenvalues for the gapped graphene and complex energies for resonances were 
derived and analysed. Exact formulae for the scattering and transfer matrices were obtained. Important features of 
the electronic states in the presence of the asymmetric potential were found in ??. Our goal in this work is to study 
the impact of these features on the electronic transport in the monolayer graphene. 

II. BASIC EQUATIONS 

The Dirac equation describing electronic states in zero-gap graphene reads @ 

^-ivph cTf^dfj, - aaSmvp^ = (E -V)-ip, (1) 

where vp is the Fermi velocity, arc the Pauli matrices, (r) is the two-component spinor.The spinor structure 
takes into account the two-sublattice structure of graphene 5m (r) and V{r) are the local perturbations of the mass 
(gap) and the chemical potential. A local mass perturbation can be induced by defects in the graphene film or in the 
substrate. The perturbation matrix elements 

diag{Vi,V2)ro5{r - ro) (2) 

are related to the a, b parameters as follows 

-Vi=a + b, -V2=a^b (3) 

The delta function perturbation is the simplest solvable short-range model. Finite radius tq plays a role of the 
regulator and is necessary in order to exclude deep states of the atomic energy scale. The finite perturbation radius 
rg leads to the quasi-momentum space form- factor proportional to the Bessel function that justifies our neglect of 
transitions between the Brillouin band points K and K' . 

Let us introduce dimensionalless variables and parameters: 



hvp/ro^ hvp/ro^ hvp/ro' ' ' ro ' ^ 9r 

roS{r - ro) = S(r - 1). (5) 



The equation ([T]) takes the form 



i ^ cr^9p + (Tarn - (TsSm \ ^ {e - u)ip, (6) 
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We confine ourselves here by the case of zero-gap graphene; m 
Let us present the two-component spinor in the form 



0. 



exp {—iEt) 



fj (r) exp [i {j 
9j (r) exp [i (j 



1/2), 
1/2), 



(7) 



where j is the pseudospin quantum number; j — ±1/2, ±3/2,.... In the opposite to the relativistic theory, this 
quantum number has nothing to do with the real spin and indicates the degeneracy in the biconic Dirac point. The 
upper fj (r) and lower gj (r) components of the spinor satisfy the equations 



(8) 



(9) 



The general solution can be found solving the second-order equation obtained by excluding one of the spinor 
components from the equation set ([8]), (jH)) in the domains < r < 7'o and r > tq : 



+ 



J U - 1) 



fj = 0. 



This equation is related to the Bessel one. Its general solution in the domain < r < tq reads 



(10) 



(11) 



where (z) and N^, (z) are the Bessel and Neumann functions respectively. The constant C2 vanishes in the domain 
< r < 1 since the solution must be regular at the origin. The solution in the region r ^ 1 reads 



f,^C,H^\/2 + (6?) 



J-1/2 



where Hi"^ (z) is Hankel's function. Matching these solutions at the circumference of the circle of radius r — 
{r — 1) we obtain a scattering matrix and a characteristic equation for the resonance states. Calculating the ratio of 
the out-going and in-going waves, we obtain the S-matrix components in the angular momentum representation (for 
details see 



jr(2) 

where -^j""* is given by the formula: 

-^j"^ = (^.-1/2 (^) i?j?i/2 (^) - ^. + 1/2 (^) H]%, (.)) - 

\a - h) J,+i/2 (e) i?j:l/2 (e) + {a + b) J,.i/2 (e) ^jl/2 (e) 
a = l, 2. 

Poles of the scattering matrix (|12p are determined by the characteristic equation 



(12) 



(13) 



-^r^ (^) = 0, 



or 



(14) 



(a - h) (e) + {a + h) (e) = z [(a - b) J,+i/2 (e) 7V,+i/2 (e) + {a + b) J,_i/2 (e) (e)] = (15) 
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III. ANALYSIS OF THE CHARACTERISTIC EQUATION AND CALCULATION OF THE 

CONDUCTIVITY 

Using the relations Hn^ (z) = J„ + iNn, Hn^ = J„ — iNn, we can write S-matrix in the form: 



A, (e) - zB, (e) _ B, (e) + lA, (e) 



A,{e) + iB,ie) B, (e) - ^A, (e)' 
and, therefore, it can be presented in the standard form [? ] 

Sj (e) = exp [i2Sj (e)] , 
where the scattering phase is given by the expression 



Sj (e) = arctan 



B,{ey 



(16) 



(17) 



(18) 



Formulae (|16p . ()17p show that the scattering matrix Sj (e) is unitary on the continuum spectrum. The functions 
Aj (e) and Bj (e) are determined as follows 



A, (e) = - {a + b) (e) + {a-h) jf^,/, (e) 



(19) 



Bj (e) = - — [(a + b) J,_i/2 (e) iV,_i/2 (e) + {a ~ b) J,+i/2 (e) (e)] + 



ne 



[jj+i/2 (e) iV,-i/2 (e) - Jj-i/2 (e) (e)] 



(20) 



Asymptotic behaviour of the scattering_phases and other scattering data at e — ^ can be obtained expanding the 
cylinder functions for small arguments 



J (z)'^—(z/2? N(z]-l -^(2/^)" forn>0, 

Jn (z) [z/Z) , iV„ [z) - I ^2^^^ ^^^^^2) for n = 0, ^ ^ 0, 



(21) 



where 7^; « 0, 577 is the Eyler-Masceroni constant, F (n) is the gamma-function. Then we have for the scattering 
phases in the lower order in e: 



tan(5±i/2 w (fo± a) -e, e ^- 0, 



(22) 



tan <5±(„+i/2) ~ ±7r (e/2)^"+^ (5 ± a) e ^ 



(23) 



The transport cross section can be written in terms of the scattering phases [3] (we have returned to dimensional 
variables here) 



2ro 



^^^^ ("^i+i ~ '^i) • 



(24) 



The transport relaxation time can be calculated using the following relation: 

l/rtr = NiVpEtr. (25) 

Taking into account the approximate formulae for phases (|22|) . ([23]) . the series ([24]) can be written in the following 
asymptotic form for e — > 0: 



2ro 



{Si/2 - + {S3/2 - + (^-3/2 - S-1/2)'^ 



(<5„+i/2 - '^n-1/2)^ + {S-n-1/2 - <5_„+i/2)^ + • • • ~ eroTT^ \2a^ + O (e)] . 



(26) 
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Then the asymptotic formula for the transport relaxation time reads: 

XjTtr tNiVF-K^la^ [1 + (e)] (27) 

It is seen from (j27p that asymptotic behaviour of the relaxation time at e —> is determined by the parameter a, i. 
e. by the symmetric component of the perturbation. 

Let us consider now the Born approximation for the scattering amplitude. The partial wave series for the transport 
cross section converges rather slowly. That is why we consider a behaviour of the transport cross section without 
the partial wave expansion. In return we can use the Born approximation in this limit. The Born formula for the 
scattering amplitude reads Ifl]: 

/-™fe^) = -^/^^(q), (28) 
where flq = /zp — Vc^ is the transferred momentum, q — 2psm^/2, V (q) is the perturbation Fourier transform: 

d^re-'^^^V {r) ^ drrV {r) d<j)exp[-iqr cos (j)] ^ 2n drrV {r) Jq {qr) . (29) 
Jo Jo Jo 

Inserting the potental 1^ into ([^^ we obtain 

V, (q) = V (p, 0) = 27rK°r2 Jq (2pro sin 0/2). (30) 
Substituting ([50)1 into we obtain the scattering amplitude: 

/f °™ {p,0) = (2prosin0/2) , (31) 

Now we can calculate the transport cross section Q: 



Sf,°™=/ c^0(l-cos0)|/^°™(p,^^)| ={pro)ro[^-^] 7r/2 / d0 (1 - cos0) (2prosin0/2) . (32) 
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lo yfivp/ro. 
This integral can be expressed in terms of the hypergeometric functions Q : 

^ d0{l- cos 0) 4 (2pro sin 6/2) = T ^ ^ • 

3F4 (3/2, 1/2, 1, 2, 1, 1, 1; - (2pro)') , (33) 

whereT^^I' ^ = r(Ji)"r\fe^)r(fe) ' ^ (o^) gamma function, 3F4 {ai,a2,a3; Pi, (32, fis, Piix) is the gener- 

alized hypergeometric function. It is determined by the series 0: 



(«l)fc («2)fc («3)fc x'' 

e,(/3i)fc(/?2),(/33),(/34)fc A:!' 



3F4 (ai, a2, as; /?!, /32, /Ss, l^r, x)^^ y^VT^T^^^V?^ TT ' (34) 

fc=0 



where (a)j, = is the rising Pohhammer symbol. When pr^ < 1, we can neglect all higher terms of this series 

and otain in result the transport scattering cross section 'E.f°™~prQ. Notice that the Born approximation is here 
asymptotically exact in the limit pr^ 0. It must be noted here that the limit of pro, Ep, ksT, fko tending to zero 
is obviously nontrivial and many-particle effects must be taken into account in this case. This problem will be not 
discussed here. The opposite limit of pro > 1 is not actual for scattering on rare point defects. It can be urgent for 
the model of random potential with large correlation radius and for the quantum dot model, but they are outside 
of the scope of this work. However, there exist at least two other possibilities to obtan the cross section behaviour, 
different from the power law 'E.f°'^^~prQ. First of them corresponds to a random potential with the correlation radius 
rg >> a. In this case one can consider scattering of electrons with prg ^ 1. Another possibility is related to the 
resonance scattering studied in our works We analyse here the latter case. 
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Let us consider the limit of large angular momentum j » 1 using the Bessel function asymptotics: 

J.(z)^-^(g)'^(l + 0(lH), ^^oo, (35) 



(36) 



(1 + 0(1/:/)), ly-^oo, 



(37) 



where e is the base of natural logarithms. The characteristic equation [15] takes the following form in the limit of 



a — b 



(2j + 1) (2j - 1 



,2j+2 



a + 6 



(2j ~ 1) (2j - 1) 



2j 



a — b 



a + b 



2j-1Vj + 1/2 J -1/2; 2j-l 



(38) 



Applying the same approach to the scattering 

Our numerical analysis shows that the energy dependence of the relaxation time crosses over to an approximate 
constant at high energy, and position of the crossing-over point depends on the ratio a/b (see fig. 1). The Boltzmannian 
conductivity is determined by the formula: 



where Ttr is determined as follows 



The mobility can be determined as the ratio: 



cr = — {EpTtr/K) , 



(39) 



(40) 



a 

en' 



where the carrier density at low temperature is determined as follows 



N/S 



2lT \hvF 



(41) 



(42) 



Using (P?]). (15^ . and PO]) we conclude that the conductivity tends to a constant value in the limit Ep ~^ O.This 
limit must be considered taking into account the fact tha the accepted here Boltzmann kinetics in invalid in the 
vicinity of the point Ep — 0. In the opposite limit of large Ep, we conclude using (??) and (^0]) instead of (l?fl) that 
conductivity increases linearly with the Fermi energy in this energy region. The mobility behaves at small and large 
Fermi energy respectively as Ep'^ and Ep^. 
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